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Bootstrap Approximations in Contractor Renormalization 
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We propose a bootstrap method for approximating the long-range terms in the Contractor Renor- 
malization (CORE) method. The idea is tested on the 2-D Heisenberg antiferromagnet and the 
frustrated J2 — Ji model. We obtain renormalization group flows that directly reveal the Neel phase 
of the unfrustrated HAF and the existence of a phase transition in the J2 — Ji model for weak 
frustration. However, we find that this bootstrap method is dependent on blocking and truncation 
schemes. For this reason, we discuss these dependencies and unresolved issues that researchers who 
use this approach must consider. 

PACS numbers: 75.40.Mg, 75.50.Ee, 02.70.-c 



I. INTRODUCTION 

Contractor Renormalization (CORE) [ij is a numerical 
method for finding the low energy spectrum and generat- 
ing the renormalization group flow of operators in quan- 
tum many-body models. Although it has been applied 
successfully to many systems (see references in Ref.Q, 
and for recent works, Refs.jl, Q), it relies on a cluster 
expansion whose higher order terms are, in most cases, 
difficult to compute. However, as we noted in Ref. in 
order to have confidence that a CORE computation is 
reliable, the convergence of the cluster expansion has to 
be checked carefully. This is particularly true for higher- 
dimensional systems. Unfortunately, for even the sim- 
plest 2-D systems, the computation of operators beyond 
nearest-neighbor blocks is a formidable task since, if we 
suppose the number of states per block is m, and the 
cluster expansion up to the fc-th term is defined to in- 
clude all connected operators in a k x k configuration of 
blocks, then the brute force calculation of the fc-th term 
would require exact diagonalization of an m*^ x m'^ ma- 
trix. Clearly this is not feasible for large values of fc. This 
paper reports our attempts to devise schemes for approxi- 
mating these longer-range terms in the cluster expansion. 

Getting a reasonably reliable handle on longer-range 
terms is important for two reasons. First, we would like 
to obtain the best accuracy for the ground-state and ex- 
cited state energies, as well as the best values for criti- 
cal exponents and other order parameters. Second, and 
possibly more important, higher order terms will gener- 
ate long-range operators that can introduce frustration, 
and therefore possibly produce qualitative changes in the 
renormalization group flows. In this paper we propose a 



CORE bootstrap scheme, wherein we use a variation of 
CORE itself to compute approximations to these long- 
range terms. The key idea behind this approach is that, 
since our focus is on the computation of new connected 
operators for a small number of blocks, we use a block 
adapted version of CORE to approximate the cluster 
Hamiltonian. By block adapted we mean that we use 
a truncation and renormalization procedure specifically 
tailored to the finite size and geometry of the particular 
configuration. 

The outline of this paper is as follows. In Section [Til 
we first test a variety of bootstrap schemes for the case of 
an unfrustrated antiferromagnet (HAF), since the ground 
state energy density of this model is known exactly [![. 
We show that the increase in accuracy we obtain for 
the ground state energy density can show significant de- 
pendence upon the blocking scheme. Fortunately, all 
schemes consistently reproduce a renormalization fiow 
that illuminates the Neel nature of the HAF ground state. 
We should remark that these schemes can be used to cal- 
culate observables such as the staggered magnetization in 
the manner shown in Ref.j2l- Adding observables to the 
Hamiltonian, however, typically breaks some symmetries 
and significantly complicates the computation. There- 
fore we only focus on the energies and the couplings in 
this exploratory work of the bootstrap idea. 
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In Section Hill we take the best blocking scheme from 
the previous section and apply it to a specific frustrated 
antiferromagnet, the J1 — J2 model. This scheme predicts 
a renormalization flow as a function of frustration, in 
which a phase transition appears around 0.37 < J2/ Ji < 
0.38. This is in rough agreement with literature, and it is 
the first time such a renormalization picture is obtained 
for this well-studied model. Section |IV] reports on the 
results of some variations of this bootstrap procedure and 
the puzzles they pose. Finally, we conclude Section |V] 
with a discussion of what we believe the implications of 
these results to be. 
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II. BLOCKING SCHEMES FOR THE 
BOOTSTRAP STEP 

As in Ref. [2] , we begin our study of the CORE boot- 
strap by studying the 2-D spin-1/2 Heisenberg Antifer- 
romagnet: 



H 



(1) 



<jj> 



As before, we assume nine-site square blocks for the basic 
CORE renormahzation group computation and truncate 
to the two lowest lying eigenstates of the corresponding 
nine-site Hamiltonian. The simplicity of the lowest range 
CORE computation for this system makes it a good test 
bed for a computation of the four-block plaquette oper- 
ator. Note that the CORE prescription for computing 
this operator exactly would be require computing those 
lowest lying eigenstates of the 2'^^ x 2^^ cluster Hamil- 
tonian which have a non-vanishing overlap with the 16 
tensor product states constructed in the truncation step. 
In general, this means computing more than the 16 low- 
est lying eigenstates of the 36-site problem. This is a 
formidable problem and it cannot be done on a PC. 

A simplification which follows from the fact that our 
renormahzation procedure preserves the rotational in- 
variance of the original theory, is that the renormalized 
Hamiltonian must be rotationally invariant. This means 
that, since we truncate to the lowest spin-i state per 
block, the generic two-block and the generic plaquette 
operators that appear in the Hamiltonian will have the 
form: 

^2 = Cuu + ^C^xSl ■ S2 (2) 
hi — Cunuu^ + 'iCxxuu{Sl ■ S2 + S2 ■ S3 

+S3-S4 + Si-Si)+4c xuxu{Sl ■ S3 -h S2 • >5'4) 

+ l6Cxxzz{Sl -82 S3 ■ S4 + Si ■ S4 S2 ■ S3) 

+ i6cxzxziSi -53 52 • 54) (3) 

where the 36 sites in the four nine-site blocks are labeled 
as shown in FigU] The constants appearing in Eq. [3] 
are put in to ensure that the coefficients c^zxz are in 
front of operators whose trace norm is unity. The coef- 
ficients Cxuxu, Cxzxz indicate the strength of diagonal in- 
teractions, while Cxxum Cxxzz are for horizontal/vertical 
interactions. Note that both the Cxx and Cxxuu terms 
are S ■ S on two neighboring sites, so the total horizon- 
tal or vertical coupling should be {c^x + '^Cxxuu)- Since 
the overall scale is irrelevant when we study the flow of 
couplings, in the plots to follow we will often express the 
couphngs in Eq. [3] in units of {cxx + '2-Cxxuu)- Given this 
notation, we now turn to a discussion of several possible 
bootstrap schemes: the Whole Block Buffering scheme, 
the Three-site Strip Blocking scheme, the Contour Line 
Blocking scheme and the "8+1" Blocking scheme. 




FIG. 1: The labeling of the four nine-site blocks. 



A. A First Attempt: Whole Block Buffering 

The key to the CORE bootstrap is the introduction of 
an intermediate, or buffering, step between each of the 
original CORE renormahzation group steps. In this in- 
termediate step we retain, along with the original spin-^ 
doublet retained in the first step, additional spin multi- 
plets so that the final number of retained states per block 
is 2^ > fc > 2. This means that for a 4-block cluster there 
will be a total of /c^ retained states. 

Next, we use the CORE contraction procedure to con- 
struct a fc** X fc"* Hamiltonian which approximates the 
full 2'^^ X in the following sense: one its low lying 
states should have nearly the same eigenvalue spectrum 
as in the full 36-site problem; two, if one restricts to the 
lowest two states of the single block Hamiltonian, then 
the renormahzation group step should produce a reason- 
able approximation to the connected operators we would 
obtain by doing the exact 36-site computation. To be 
precise, this approximating fc^ x Hamiltonian is con- 
structed as follows: first, we compute the k x k range-1 
(single block) terms and k^ x k^ connected range-2 (two 
adjacent block) terms; second, we add these terms up 
to construct the approximate Hamiltonian. The hope is 
that since, if we kept all 2^ states per block, there would 
be no plaquette operator at all, keeping the range-2 ap- 
proximation for fc > 2 states will, for big enough k, be 
a reasonable thing to do. The only question is how big 
does k have to be for this to be true. 

An important feature of this procedure is that we 
keep the original spin-^ multiplet as one of our k re- 
tained states. This means that, if we carry out a nearest 
neighbor renormahzation group step for this approximate 
Hamiltonian, keeping only the lowest spin-i multiplet 
for a single site, then we will automatically reproduce 
the result obtained in the full eighteen-site, two-block 
calculation. Obviously, this procedure can be general- 
ized in ways that fail to guarantee that the connected 
range-1 and connected range-2 terms computed for the 
effective problem will be the same as those obtained from 
the original 18-site calculations. In that case, however, 
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FIG. 2: Approximate energy density as a function of k, the 
number of states per block in the buffering step. The four k 
values are chosen to preserve whole spin multiplets in the first 
iteration. For subsequent iterations, we keep instead all the 
whole multiplets that can be fit within the k lowest states. 



the question of what is the proper definition of the con- 
nected plaquettc operator raises its ugly head. To avoid 
this complication, in what follows, we will only discuss 
buffering stages for which this cannot happen. 

The purpose of the buffer step is to reduce numeri- 
cal effort, so we limit ourselves to calculating range-2 
configurations - those connecting horizontal, vertical and 
diagonal nine-site blocks. A point worth noting is that 
after the first CORE iteration our renormalized Hamilto- 
nian will contain a connected plaquette term which cou- 
ples every four site square. Limiting ourselves to range-2 
means the c^xzz , Cxzxz terms at the very center of the 
36-site cannot be included. While these operators have 
very small norms, this approach also violates the diame- 
ter expansion rule proposed in Ref. by not including all 
diameter- \/2 operators. This is potentially a defect of the 
whole block buffering scheme. Nevertheless, for the pur- 
pose of testing and comparison, we follow this strategy 
and see what happens for several values of k. The re- 
sulting energy densities are shown in FigO The general 
question of how best to perform whole block buffering 
without discarding any terms is an issue that we may 
return to at a later date. 

Fig[2] shows that as we add more states, we get closer 
and closer to the exact energy density. This is nontrivial 
and appears to indicate that we are on the right track, 
but notice that for fc = 6, the energy density at '-^ —0.68 is 
much worse than the —0.666 obtained without bootstrap 
. This may be the result of the defect mentioned above. 

Apart from the energy density, this calculation (as well 
as all other blocking schemes to be discussed) reveals an 
important feature of the HAF. With each successive it- 
eration, the sum of Cxx + '^Cxuxu goes to zero and Cxuxu 




FIG. 3: Two decoupled inter-penetrating lattice. The fact 
that the HAF runs toward this configuration illustrates the 
Neel character of the ground state. 
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FIG. 4: Three-site strip blocking. Shown on the right are 
operators for the buffering step that can be embedded in the 
36 sites that involve no more than the diagonalization of two 
original nine-site blocks. 



becomes more and more negative. At the third iteration 
of fc = 14, for example, we have the following coefficients: 

C-xuxuli^^xx ~^ ^^xuxu) ~ ''^4, Cxxzz/i^^xx ~t- '^Cxuxu) ~ 

0.006, Cxzxz/ [cxx + '^Cxuxu) ~ -4. This means that un- 
der the renormalization group flow our lattice Hamilto- 
nian flows to that of two weakly-coupled interpenetrating 
ferromagnets (Fig[3]). This picture appears to provide a 
very attractive qualitative explanation as to why, at long 
range, the HAF ground state is Neel-like. 



B. Three-site Strips 

The previous blocking scheme leaves the nine-site 
blocks intact, but nothing prevents us from breaking 
them up in the buffering step. Appropriately dividing 
the block into smaller pieces allows us to include more 
interactions at a smaller computational cost. More im- 
portantly, this flexibility allows us to retain four-body 
interactions we previously discarded, and so we will do 
not have to violate the idea of ordering the cluster ex- 
pansion by diameter/proximity. 

Fig[3] shows a way of breaking the nine-site block into 
three rectangular strips, along with a list of operators 
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FIG. 5: (Color Online) Running of horizontal/ vertical cou- 
plings in four CORE iterations using alternating strip block- 
ing. The two-body horizonal coupling is defined to be 
2 (c^ +2cS^„„ ) / (cL +2cS^„„ -l-c^^ -|-2c^^„„ ) and the four-body 
horizontal coupling is 2c^^^^/{c*l^ -|- 2cS^„„ -I- cl^ + 2c^^„„). 
The vertical couplings are defined similarly. 



calculated in the buffering step. Embedding these terms 
in the 36-site problem, we see that there are two six-strip 
terms that connect all four nine-site blocks. For this strip 
blocking we keep two states (a spin-i) per strip, which 
means that a six-strip term requires contracting the low- 
est of 2^* states to 2^ states. Computationally this is 
not much more demanding than the 2^^-to-2^ contrac- 
tion required for the original range-2 calculation without 
bootstrap. Thus, in this way we can take into account 
the four-block interaction more satisfactorily than in the 
whole block buffering scheme. 

The obvious problem with this approach is the loss of 
rotational symmetry. Now, the up-down couplings will, 
in general, be different from the left-right couplings. This 
difference can and will grow as we run the renormaliza- 
tion group flow. This, of course, limits the number of 
RG steps which can be carried out and so, if the limiting 
behavior of the renormalized Hamiltonian hasn't clearly 
emerged before these asymmetries grow too large one 
loses the ability to understand the long distance behavior 
of the theory. We can mitigate this problem somewhat by 
alternating the blocking direction at every iteration, but 
this only improves the energy from —0.6678 to —0.6679. 
Compared to the exact value at —0.6694. Both numbers 
are ~ 0.2% accurate, as was the fc = 14 result in Fig. [21 

FigO shows how the asymmetry between horizontal 
and vertical couplings grows with renormalization un- 
der the alternating strip-blocking scheme. Just as in the 
whole-block case, both horizontal and vertical couplings 
are eventually dominated by the diagonal couplings. To 
see this more clearly, we artificially restore the rotation 
symmetry by averaging the horizontal and vertical in- 




Number of iterations 



FIG. 6: (Color Online) Running of diagonal couplings in 
four CORE iterations using strip blocking with averaging of 
coefficients. The two-body diagonal term is 2c^^^^/{c^x + 
'2Cxxuu + + 2Cxxuu) and the four-body diagonal term is 
'^Cxzxz/ i^^xx "t" '^Cxxuu + *^xx "t" '^^xxuu)- This shows that the 
HAF runs to the two-lattice limit in Fig[3] 

teraction, i.e. set cj^. ^ (cj^. + d^x)/'^ ^tc. where the 
superscripts indicate horizontal and vertical coefficients. 
The averaging scheme hardly changes the energy density, 
which now becomes —0.6677, but here we get a clean pic- 
ture (Figini) of the growth of diagonal couplings. The flow 
to two weakly-coupled interpenetrating ferromagnets is 
once again made manifest. 



C. Contour Lines 

Since three-site strips artiflcially break the rotational 
symmetry, one may wonder if we can rearrange the strips 
to avoid this. Fig. [7] exhibits two arrangements that 
achieve this by breaking up the nine sites into 1-1-3 + 5 
sites, forming a contour line structure. Notice that this 
is no longer a scheme which could be used to compute 
an RG flow for the full lattice, because no matter how 
we orient the contour lines, no single arrangement covers 
the entire lattice (four arrangements would be needed). 
Nevertheless, in the buffering step there is no barrier to 
using either of these arrangements to construct the ap- 
proximation to the 36-Hamiltonian. 

Table U lists the results of a number of different contour 
line blockings. Blocking scheme D) makes use of the 
racetrack contour in Fig|7]and keeps a spin-^ multiplet 
in each strip (made of one, three or flve sites), which is 
always possible with an odd number of sites. This scheme 
does not include any extra long-range terms apart from 
the term connecting two adjacent nine-site blocks (made 
of six strips), so it is similar to the whole block buffering 
scheme. However, instead of including the interaction 
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FIG. 7: Two different types of rotationally symmetric con- 
tour line blocking. We will refer to the left configuration as 
a racetrack contour, and the right as a cross contour. For 
the racetrack contour, we have calculated two operators that 
couples all four nine-site blocks; the 4x4 square in the mid- 
dle (shaded above), and the 20-site outer ring. For the cross 
contour, the 20-site cross in the middle is the term connecting 
four blocks we can handle. 



TABLE I: Performance of strip and contour-line blocking as 
measured by energy density. Two states per strip is assumed 
except for G) 



Bootstrap Scheme Used 


Approx. Energy Density 
and Percentage Error 
(Exact = -0.66943) 


A) Three-site strip 


-0.6678 (+0.24%) 


B) A) with alternating 
orientation 


-0.6679 (+0.22%) 


C) A) with horizontal/vertical 
coefficients averaged 


-0.6677 (+0.25%) 


D) Racetrack contour without 
extra long-range terms 


-0.6729 (-0.52%) 


E) D) with 4x4 square in the 
middle included 


-0.6658 (+0.54% 


F) E) with 20-site outer ring 


-0.6646 (+0.72%) 


G) Like D) but with four 

states instead of two kept in 
each five-site strip 


-0.6669 (+0.37%) 


H) Cross contour with 20-site 
cross in the middle 


-0.6718 (-0.36%) 


I) Coefficients from F) 

and H) averaged at the end 
of each iteration 


-0.6668 (+0.39%) 



between two diagonal nine-site blocks, we only include 
the four-site interaction in the middle, and this does not 
have to be calculated because the four sites in the middle 
have been kept free. 

Blocking scheme E) takes this further by adding the 
operator corresponding to the 16 sites in the middle. To 
make sure we are only adding connected operators, we 
have to calculate embedded sub-blocks, such as the 2x4 
rectangle, etc. Blocking scheme F) includes yet another 
term corresponding to the 20-site outer ring. Note that 
there are still many long-range operators that can be em- 
bedded in the 36 sites under racetrack contour that we 
have not included. For instance, we have not calculated 
the 20-site "T" -shaped term corresponding to two nine- 
site blocks plus two other free-sites in the middle. How- 
ever, it would seem that E) and F) already encodes most 
of the interaction coupling the four nine-site blocks, and 
this we know is important because the diagonal coupling 
grows with renormalization. 

As apparent from the energy density, our hope was not 
borne out. The energy densities obtained from E) and 
F) are even further from the exact value than the result 
without bootstrap To see whether including more 
states would help, we modify D) by increasing the num- 
ber of states in five-site strips from two to four (i.e. two 
spin-i multiplets). This is blocking scheme G), and while 
it improves the accuracy, it is only about as good as the 
non-bootstrap calculation. Switching from the racetrack 
contour to the cross contour and including the 20-site for 
the strips in the middle (shaded on the right of Figl?]) 
gets us about the same level of accuracy. Finally, artifi- 
cially averaging the coefficients from the racetrack con- 
tours with those from the cross contours does not help 
either. 

This experiment appears to teach us two things. First, 
we see that the numbers roughly bracket the exact value, 
so just as in Fig[21 this indicates that a good approxi- 
mation of the four-block plaquette should yield a very 
good energy density. On the other hand we learn that 
a complex blocking scheme is prone to over-emphasizing 
certain interactions, and that adding longer-range opera- 
tors in such situation does not necessarily lead to better 
accuracy. In absence of something resembling the diam- 
eter expansion rule Q (which only makes sense in large, 
uniform lattices), it is not clear what operators we should 
calculate to get improvements. 



D. The "8+1" Blocking Scheme 

With the lessons of our previous attempts in mind, 
we now present a blocking scheme that produces much 
better accuracy. We saw that the whole blocks in Section 
III Al gave reasonably good results and did not assume 
arbitrary structures. The problem there was that the 
four-site operator in the middle, which would be present 
after the first iteration or in a frustrated model, cannot 
be handled in a simple, natural manner. The minimal 
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FIG. 8: The "8+1" blocking scheme. One site out ol each 
nine-site block is left free, enabling the interaction on the 
four sites in the middle to be included naturally. 

way to deal with this is to break the nine sites into 8+1 
sites, leaving the corner site towards the middle of the 
four blocks free. This is shown in FiglH 

To mimic the 14-state whole block calculation, we keep 
all the multiplets that can be fit within the lowest eight 
states in the L-shaped eight sites. This turns out to be 
seven states in the first iteration, so together with the free 
site we would have 14 states per nine-site block. This 
simple scheme yields a remarkable value, —0.66927 as 
energy density, which is +0.03% from the exact value. 
With this encouraging result, we are ready to see what 
happens in the more complicated J\ — Ji model. 



III. THE Ji Ja MODEL 

As we have seen, CORE's way of exhibiting the long- 
range Neel order of the HAF, is that RG with bootstrap 
always generates strong diagonal couplings; i.e., in all the 
blocking schemes Cxuxul{c-xx + ^^c^^xun) flows to a large 
negative number. This raises the question of what hap- 
pens if we add a large positive c^uxu at the initial step. 
This, of course, brings us to the well known frustrated 
Ji — Ji model: 

H ^ J^^S -S^ ,h^S -S. (4) 

nn diag 

where the first sum is over nearest-neighbors and the 
second sum is over diagonal sites. This corresponds to 

\^xx — Jl:Cxuxu — J2lJl^^xxuu — C-xzxz — C-xxzz — 

0}. The J2 — Ji model has been extensively studied in 
the condensed-matter literature ([1, Hill and references 
therein) and it is generally believed that there are sev- 
eral phase transitions within the region < J2/J1 < 1. 
Ref. [6| for example lists four conjectured critical points 
at J2IJ1 = 0.34, 0.38, 0.5, 0.62. In this work the ground 
state transitions from the Neel phase to various dimerized 
configurations with increasing frustration . 




FIG. 9: (Color Online) Renormalization group flow generated 
by "8+1" bootstrap CORE for weak frustration J2/J1 = 
{0.2, 0.3, 0.35, 0.371, 0.372, 0.373, 0.374, 0.375, 0.39, 0.4, 0.45}. 
All the couplings above are implicitly divided by 
{cxx + 2cxxuu), SO that the overall scale does not affect 
our analysis of phases. We run the system starting from 
several values of Cxuxu from 0.2 to 0.45; each arrow on the 
lines indicate a CORE iteration, their size indicates amount 
of change in each step. To avoid cluttering, this region in 
parameter space has been selected to give a clear picture of 
the quantum phase transition. 



Past works on this subject typically rely on two types 
of techniques. One is to use perturbative methods around 
some conjectured starting configurations; the other is to 
perform calculations on finite lattices with some numer- 
ical methods and extrapolate to an infinite lattice. The 
bootstrap CORE is therefore a very different method and 
may provide independent checks on past results. Since 
we have an RG flow picture that describes the Neel phase 
at J2/ Ji = 0, we would expect to find at least one critical 
point where this picture begins to change. 

In what follows, using the "8+1" bootstrap scheme, 
we compute the renormalization flow for J2I J\ < 0.5, as 
shown in Fig[51 For J2/J1 < 0.371, the system flows to- 
wards the unfrustrated case, whereas J2/J1 > 0.372, the 
frustration grows. The qualitative feature of the entire 
region can be seen quite clearly from the plot. To our 
knowledge this is the first time such a renormalization 
group result has been obtained directly from first prin- 
ciples for the J2 — Ji model. In Ref.ff?], J2/J1 = 0.34 
is thought to be the point where columnar spin dimers 
begin to appear. Another critical point at J2/J1 — 0.38, 
which is more established, is thought to be where the sys- 
tem becomes a disordered spin-liquid. In between these 
two values is a region where columnar dimer order and 
Neel order coexist. The critical point we see in Fig 15] 
appears to be in agreement with literature. 

Things do not go as well for J2/J1 > 0.5. As the 
frustration increases, there is more competition among 
states and the gaps decrease and the gap between the 
lowest spin-i and other excited states in the nine-site 
block become very small. In fact, after one iteration with 



0.2 





FIG. 10: (Color Online) Renormalization group flow starting 
with strong frustrations. The gap between the lowest spin-^ 
and other states in the nine-site block is often small or even 
negative in much of this region, so our algorithm may not be 
reliable here. No clear pattern emerges compared to Fig[9] 
and there are many large erratic jumps from strong to weak 
frustration region. 



J2/J1 = 0.59 or J2/J1 — 0.75, we arrive at a system 
where the nine-site block no longer has a spin-i as its 
ground state, and so the truncation scheme we are using 
should change. If we don't do this then we see that, 
starting with J2/J1 > 0.372, the frustration eventually 
grows to a point where the CORE recipe breaks down. 
FigfTUl shows this happening to the flow of couplings at 
J2/J1 > 0.5. 

While the first iteration always takes us to higher 
frustration for 0.372 < J2/J1 <~ 0.9, the flow is 
quite erratic in the strong frustration region. This pre- 
sumably occurs because we fail to modify our trun- 
cation procedure as needed. Another curious behav- 
ior that is probably also due to our failure to follow 
the true ground state, is that there is a neighborhood 
around Cxuxu 1 where the flow takes us back to 
weakly frustrated regions. For example, if we start 
with {cxuxu,Cxxzz,Cxzxz} = {1.085,-0.064,-0.223}, we 
would be taken to {0.318, -0.021, -0.086} at the next it- 
eration, and small changes at the beginning lead to very 
different couplings later on. To give an idea how seri- 
ous the gap issue is, we plot in Fig[Tl] the size of the 
first excited state gap in the nine-site block as a func- 
tion of Cxuxu/ {^xx ~^ '^Cxxuu) ^ Cxzxz / (,Cxx ~t~ ^Cxxuu') with 

c-xx = 1, Cxxuu = —0.2, Cxxzz = 0. All of these results are 
consistent with what has been seen in one-dimensional 
examples when one fails to modify the space of retained 
states as needed whenever level crossings occur. 

Compared to the couplings, we see that, because it is 
largely determined by the first two iterations, the energy 



FIG. 11: (Color Online) Plot of the gap between the low- 
est spin-i multiplet and the first excited states in the nine- 
site block in a typical frustrated region. The couplings are 
implicitly in units of {c^x + 'ic^xnu), with c^x = '^,Cxxuu = 
—0.2, Cxxzz = 0. The gap is set to zero when the spin-i mul- 
tiplet ceases to be the ground states. 



density is a little less sensitive to the small gaps. An 
plot of the energy density with respect to J2 / Ji is shown 
in Fig[T21 This can be contrasted with Fig. 5 of Ref.jl], 
where a qualitatively similar plot is obtained from finite 
lattice extrapolation. Instead of a peak and discontinuity 
at J2/J1 ^ 0.6, we find a peak around 0.53. No partic- 
ular structure is found near J2I J\ ^ 0.6, but the slope 
appears to be the steepest there. It is possible that these 
features are indirect signs of conjectured phase transi- 
tions at J2/J1 = 0.5,0.6 (similar features are found in 
finite lattice entropy plots [^), but we would have to 
run magnetization and other specific operators such as 
columnar dimers [6] to characterize the phases in detail. 



IV. PERFORMANCE OF OTHER VARIANTS 

In conclusion, we would say that our bootstrap algo- 
rithm appears to be a promising way of increasing the ac- 
curacy of CORE computations. However, we don't have 
a general prescription for estimating the number of states 
to keep in the buffering step. Moreover, we haven't yet 
tried to do the more difficult questions of how to handle 
the truncation scheme when we flow to points where level 
crossings occur. 

The issue of strong frustration aside, there are other 
questions which need more study. How robust are our 
results? What speciflc features in the blocking scheme 
work best? How do we characterize the different phases 
believed to exist in the J\ — Ji model, if the flow of the 
coefficients in the Hamiltonian fail to clearly distinguish 
them? In this section we will discuss some preliminary 
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FIG. 12: Energy density as a function of frustration in the 
Ji — J2 model. The slope appears to be the steepest around 
0.6. The last data point at 0.75 should be lower than it ap- 
pears since we have to stop after one iteration. 

results on variants of the algorithm and the puzzles they 
pose. 

A. Use of Reduced Density Matrix 

Under the 8 + 1 bootstrap scheme the retained states 
on the eight sites are chosen according to their ener- 
gies. One might question the arbitrariness of favoring 
the energy spectrum of this irregular eight-site config- 
uration. One alternative is to adopt the reduced den- 
sity matrix truncation technique from the DMRG 
literature. The use of reduced density matrix in CORE 
first appeared in Ref .JlH . where it was suggested as a 
diagnostic tool. Ref. [2] explored this idea further and 
argued that while this diagnostic tool is not necessarily 
reliable (since CORE does not rely solely on preserva- 
tion of state structure - see Ref. [2] for comparison be- 
tween DMRG and CORE), it can provide an alternative 
truncation method. This method has yielded interesting 
results in Refs. so it is instructive to try it on the 

8 -f 1 bootstrap calculation. 

The truncation method, in the present context, means 
the following. We first consider a site configuration (a 
superhlock if we follow DMRG terminologies) that con- 
tains the eight sites of interest, i.e. a configuration that 
is the eight sites plus some "environment" sites. On the 
superblock, we will compute a target state, which can 
be the ground state of the superblock or a mixed state 
formed by some low energy states. Denoting the target 
state density matrix as a and the environment as E, we 
compute the reduced density matrix on the eight sites as 
follows: 

p = Tr^cr. (5) 



(If the target state is a pure state = ji/') (V'lO The 

eigenvectors corresponding to the largest eigenvalues of 
this matrix would tell us what states in the eight sites 
contribute the most to the formation of jV'), and so this 
gives us a way to choose what states to keep. 

It should be apparent that the result we get would 
strongly depend on how we choose the target state (and 
the environment). To make things as simple as possible in 
our test, we choose the superblock to be a 16-site square 
and the target state to be its spin-0 ground state. This 
way we don't have to worry about degeneracies or en- 
forcement of spin symmetry. Approximately seven eigen- 
states of p (depending on the multiplet structure) with 
the largest eigenvalues are chosen as the retained states. 
With this particular choice of target state, the 8-)-l boot- 
strap CORE for the unfrustrated HAF gives us an energy 
density of -0.66962, which is —0.03% from the exact value 
and as good as our original recipe. 

Unfortunately, this choice does not seem to perform as 
well in identifying phase transitions in the Ji — J2 model. 
In Ref. we have seen that targeting the superblock 
ground state could lead to improvement in ground state 
energy at the expense of accuracy in excited states. Thus 
we might expect the target state choice to be related to 
a tradeoff between the ground state energy and the over- 
all physical picture. Indeed, what we observe here is a 
subtly altered RG flow. For relatively weak frustration, 
say, J2/J1 = 0.3, we see an initial flow to the left of 
FigM i-e. the less frustrated region. Yet as we enter the 
far front left region (negative Cxuxu and negative Cxxzz), 
the renormalization turns around and changes direction 
toward strong frustration on the right. 

It is possible, of course, that the situation would im- 
prove if we use a mixed target state that targets some 
states above the ground state of the superblock. This 
technique could be important for not just the choice of 
states for the bootstrap step, but also for the choice of fi- 
nal retained states in the nine-site block. We have seen in 
FigUDthat the lowest spin-i may not sufficiently capture 
the low energy spectrum under strong frustration. If we 
would like to continue to work with spin-^, the reduced 
density matrix approach may be the only option. 



B. Resummation 

There is another degree of freedom in the procedure 
that affects the choice of retained states, and that is how 
we decide which part of the Hamiltonian belongs to a 
given set of sites. Suppose we want to include all the 
operators that act on two bodies A and B, and we know 
there is an operator O4 that acts on four bodies A, B, C 
and D. For any operator O2 that acts on A and B, we 
can always decompose O4 into (O4 — O2) + O2. Hence it 
should be equally valid to say that we have an operator 
O4 — O2 acting on the four bodies, and an operator O2 
that acts on A and B. But O2 can be anything, so what 
should we choose it to be? 
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In our case, the question is whether we should set 
^xx ^ '-'xx ~t~ ^^xxuu-) '^xxuu ^ 0: i.Q. consider the Cxxuu 
term as a truly two body term instead of part of a four 
body term. This would seem like a very natural thing 
to do, but it would dramatically change the cluster ex- 
pansion. Running 8+1 bootstrap CORE with this re- 
summation yields —0.66957 for the HAF energy density, 
and the growing ferromagnetic diagonal couplings come 
out as well. However, the behavior of couplings in the 
Ji — J2 model turns out to be very different. This time, 
instead of running toward strong frustration too early, 
the system runs to Neel for relatively large J2/J1. The 
critical point shifts to above 0.4. This change may not 
seem surprising if we note that applying resummation to 
the five-site blocks in Ref. ^ also obliterates the grow- 
ing spin picture there. The problem is how we should 
interpret the effects of resummation. 

Mathematically the issue of grouping terms in the 
Hamiltonian is hardly different from the issue of selecting 
a blocking scheme, which strongly affects the resulting 
qualitative picture in CORE. One possible explanation 
for the effect of resummation is that the original group- 
ing contains some structural information of the model 
in large scale, such as what counts as self-interactions 
within a group of blocks and what counts as long-range 
interactions. Since we would like to study the flow of 
parameters in the same parameter space, it may be more 
consistent to preserve this information. Under this inter- 
pretation, resummation can yield very accurate energies, 
but we may lose control on the qualitative picture as a 
result. 



V. FURTHER DISCUSSION 



We have shown that it is possible to approximate long- 
range terms in CORE using a bootstrap approach and a 
relatively small number of additional states. In particu- 
lar, we have obtained some remarkable renormalization 
flows that illuminate the phases of a 2-D antiferromagnet. 
However, an equally important lesson that can be drawn 
from our results is that CORE is sensitive to blocking 
and truncation schemes. This suggests that the method 
itself needs further study. Moreover, it is clear that one 
should be careful about drawing any strong conclusions 
about specific applications without checking that the im- 
portant results one arrives at are common to several dif- 
ferent schemes. 

For future directions, our work should be immediately 
extended in two ways. The first is to investigate the 
effect of using reduced density matrix truncation with 
different target states at both stages of the bootstrap 
approximation. This has the prospect of resolving sym- 
metry problems and allowing us to smoothly transition 
between truncation schemes when level crossings occur. 
This would be particularly important in extending our 
results to the highly frustrated case. Second, of course, 
is that we have not attempted to study how the expec- 
tation values of various relevant operators change in the 
J2 — Ji model with increasing frustration. Such informa- 
tion is crucial to shedding more light on the structure of 
the phases. 
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